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Abstract: In this paper we continue the development of a model for strong interactions at high energy, 
based on two ingredients: CGC/saturation approach and the BFKL Pomeron. In our approach, the un¬ 
known mechanism of confinement of quarks and gluons, is characterized by several numerical parameters, 
which are extracted from the experimental data. We demonstrate that the two channel model, successfully 
describes the experimental data, including both the value of the elastic slope and the energy behaviour of 
the single diffraction cross section. We show that the disagreement with experimental data of our previous 
single channel eikonal model [6], stems from the simplified approach used for the hadron structure, and is 
not related to our principal theoretical input, based on the CGC/saturation approach. 
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1 Introduction. 

In this paper we expand our approach, of analyzing soft interactions at high energy, based on two main 
ingredients: the Colour Glass Condensate(CGC)/saturation effective theory for high energy QGD (see 
Ref.[l]) and references therein), and the BFKL Pomeron that describes both soft and hard interactions 
at high energy [2, 3]. The idea that there is only one BFKL Pomeron, which is not even a pole in the 
angular momentum plane (not a Reggeon, see Ref. [4] for details), presumes that the unknown mechanism 
of confinement in QCD of quarks and gluons, is not important, and that its influence can be mimiced by 
the determination of several parameters of the GGG/saturation approach, which depend on long distance 
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physics. As an example of such a parameter, we mention the behavior of the scattering amplitude of the 
BFKL Pomeron at large impact parameters. This contradicts the hope that confinement would lead to a 
Pomeron which is a Regge pole, (see Ref. [5] and references there in). Unfortunately, due to the embryonic 
state of our understanding of confinement in QCD, we do not yet have a theoretical tool to differentiate 
between these two approaches. Hence, we concentrate our efforts on comparing the results of our approach, 
with relevant experimental data, hoping that an evaluation will allow us to check how viable our scenario 
is, and to find the specific features where our approach differs, from one based on soft Pomeron calculus. 

Our first attempt [6] shows, that we can describe the main features of the data, but we found two results 
in our description which imply a potential problem for our approach: the result for the elastic slope, is much 
smaller than the experimental measurement at the LHC energies, and the behaviour of the cross section 
for single diffraction, which displays oscillating saturation as a function of energy. The goal of this paper 
is to show that these problematic features, occur due to our oversimplified model in Ref. [6]. In this model, 
all shadowing corrections stem from two sources: the eikonal rescattering, and the interaction of the BFKL 
Pomerons taken in the CGC/saturation approach. In this paper we develop a two channel model instead 
of the eikonal approximation, which is responsible for the low mass component in diffraction production, 
and which is an essential ingredient of all models of the high energy hadron scattering, now on the market 
[7-16]. 

The paper is organized as follows. In the next section we describe the main theoretical input used in 
this paper. This section reviews the results that have been derived and employed in our previous paper 
[6]. In addition, we include a new formula for single diffraction dissociation, which is based on the closed 
form solution to the Balitsky-Kovchegov equation[17]. In section 3 we discuss the structure of our model, 
and the phenomenological parameters that have been introduced. In this section, based on theoretical 
considerations, we estimate the range of values of the parameters. The fourth section is devoted to the 
description and results of the fit. 


2 Theoretical input 


2.1 Dressed Pomeron 


In the CGC/saturation approach (see Ref. [I] for the review of this approach), the scattering amplitude of two 
dipoles at high energy, is described as the exchange of the dressed Pomeron, which can be calculated using 
the MPSI approximation [18], displayed in Fig. 1. In Ref. [19] it was shown, that the MPSI approximation 
is valid over a wide range of rapidities: 


V < 


A, 


In 


A2 

^BFKL 


( 2 . 1 ) 


From Fig. 1, it appears that the calculations can be performed in two stages: first, to find the sum of 
‘fan’ diagrams for rapidity regions (0, Y’) and (Y’ Y); and , second, to evaluate the diagram of Fig. 1-a 
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Figure 1. MPSI approximation: the simplest diagram(Fig. 1-a) and one Pomeron contribution (Fig. 1-b). C = 
d|/47r. X 1 — X 2 = r . ^{xi +X 2 ) = b . The wavy lines describe BFKL Pomerons. The blobs stand for triple Pomeron 
vertices. 


using the following sum rules for the BFKL Pomeron exchange [20, 21] (see Fig. 1-b for all notations) 

2 p 

^ Gjp{Y -0,r,R;b) = (fr'd%'Gip{Y-Y',r,r',h-h'^ Gjp (y'Y , R,b - b''^ . (2.2) 

In Eq. (2.1) Gjp denotes the Green’s function of the BFKL Pomeron. 

The first step can be accomplished by finding a solution to the Balitsky-Kovchegov (BK) non-linear 
equation [22]. The solution has different forms in the three kinematic regions. 

1. r‘^Q‘1 iY,b) <C 1, where Qs denotes the saturation scale[20, 23, 24]. The non-linear corrections are 
small, and the solution is the BFKL Pomeron, which has the following form[2, 3] 

Gjp{Y,r,R-,b) = (u;u;*) 5 yi^e^BFKLy-is^^ ( 2 . 3 ) 

with Abfkl = 2 In 2 0:5 and D = 14C(3)q:5 = 16.828 05 . 

Gjp {Y,r,R-,b) denotes the BFKL Pomeron Green function, as the QCD coupling, r and R are the 
sizes of two interacting dipoles. Y = Ins, where s = W^. W denotes the energy of the interaction, 
and b the impact parameter of the scattering amplitude for two dipoles. 


ww* = 



(2.4) 


From Eq. (2.3) it is obvious that the BFKL Pomeron is not a pole in angular momentum, but a branch 
cut, since its Y-dependence has an additional Ins term; it also does not reproduce the exponential 
decrease at large b, which follows from the general properties of analyticity and unitarity, for the 
exchange of the BFKL Pomeron [25]. 
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r2g2 ^Y,b) rv-' 1 (vicinity of the saturation scale). The scattering amplitude has the following form [26, 
27] 

A = Gip(z) = Const (r^ (Y, , (2.5) 

where Const denotes a constant, and where the critical anomalous dimension jcr, can be found from 


X (7cr) 
1 '7cr 


dX ilcr) 

d'Jcr 


and X(7) 


2'0(1) - V’(7) - V’(l-7)- 


( 2 . 6 ) 


3. r^Q‘1 > 1 ( inside the saturation domain). The breakthrough that allows us to develop phe¬ 

nomenology based on the CGC/saturation approach, was published in [28]. In this paper a simple 
approximation to the numerical solution of the BK equation was found, which is of the form 


N^^{Gip{z)) = a (1-exp(-Gjp(z))) + {1 - a) ^ ^ 

with a = 0.65. 


Using the solution of Eq. (2.7) we calculated the diagrams shown in Fig. 1-a using Eq. (2.2) (see 
Refs. [6, 28]). The result of this calculation gives the following expression for the Green function of the 
dressed Pomeron: 


-YQ,r,R,b) = 


a^< 1 — exp (—T {Y — Yq, r, R, b)) > -|- 2a(l — a) 


( 2 . 8 ) 


TiY-Yo,r,R,b) 

1 + T{Y-Yo,r,R,b) 


+ (1 — a)^ S 1 “ exp 


1 


T{Y -Yo,r,R,b)J T (Y - Yo,r, R,b) \ ' T (Y - Yo,r, R,b) 


T O, 


where T (x) is the incomplete Euler gamma function (see 8.35 of Ref. [29]). The function T (Y — Yq, r, R, b) 
can be found from Eq. (2.2), and has the form: 

T{Y-Yo,r,R,b) = ^ Gjp (z ^ 0) = 0o (R, U, 6))= <PqS {b) , (2.9) 

where we used two inputs: r = R and = (l/(m^R^)) S {m ,b) exp (AT). The parameter A, in leading 
order of perturbative QCD, is given by A = asxilcr) /(I ~ 7cr) ■ The parameter m and the function 
S {m , b) originate from non-perturbative QCD contributions, and are parmeterized as: 

2 p 

S{mb) = where / d%S (b) = 1 (2.10) 

<pQ can be calculated from the initial conditions using Eq. (2.3). Since, we do not know these conditions, 
we will consider (pQ as an additional phenomenological parameter. 
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2.2 Single diffraction 

The equation for the single diffraction production was proposed more than a decade ago [17]. In Ref. [17] it 
was shown that the equation has the same form as the BK equation[22, 42] for the function 

G{Y,Yo,r,b) = 2N{Y,r,b) - {Y,Yo,r,b), (2.11) 

where 

adiff{Y,Yo,r) = j cfb {Y,Yo,r,b) (2.12) 

is the cross section for diffraction production, with the rapidity gap (Ygap = Y — Ym) larger than Yq{ 
Ygap > To)- 

N {Y,r,b) is the imaginary part of the elastic amplitude, other notations are clarified in Fig. 2. 

The difference between the elastic amplitude and G {Y,YQ,r,b), is only in the initial condition, which 
for G (y, Iq, r, b) is given by the following equation: 

G{Y,Yo = Y,r,b) = 2N{Y,r,b) - {Y,r,b) (2.13) 


where the last term denotes the elastic cross section. 

Bearing in mind that the solution to Eq. (2.13) is given by Eq. (2.7), we can obtain a solution for single 
diffractive production. The cross section for the production of a bunch of hadrons with a mass from Mmin 
to Mmax can be written as 


^dif f {Y, Ymax ,Ymin,r) = j dH {Y, Yrr, 


.Yr 


max^ ^ min’i ' ? ' 


(2.14) 


where the amplitude , {Y, Y^ax, Yminf] b) takes the form: 

NsD{Y,Ymax,Ymin,Yb) = {TT {Y,YmaxM) " {TT {Y,YminM) (2-15) 

where Y^ax = In and = In {mY^/sq). 

For TT (Y, Ym , b) we have the following expression 

TT{Y,YM,b) = ^2G‘"'"=="‘'(r(y-yM,6)) - (G‘'”=="‘'(r(y-yM,6)))^^ (2.16) 

2.3 Double diffraction 

We use s-channel unitarity, to obtain the expression for the cross section for double diffractive production. 
Indeed, the unitarity constraints for the dressed Pomeron given by the diagrams of Fig. 1-a take the form 


2(^dressed (2 T (E; 6)) + NDD{Y-,b) 


(2.17) 
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Figure 2. MPSI approximation: the simplest diagrams for single diffraction production. The wavy lines describe 
BFKL Pomerons. The blobs stand for triple Pomeron vertices. The dotted line denotes the cut Pomeron. Ym = 
In (M^/sq), where M is the mass of produced particles and sq is the scale taken to be of the order of 1 GeV^. 

where G'dressed ^ 27 " (Y]b)) describes all inelastic processes that are generated by the dressed Pomeron ex¬ 
change. One can check this formula for the exchange of a single BFKL Pomeron. The general proof that 
the inelastic cross section due to the exchange of the dressed Pomeron, is given by (^dressed {2T {Y]b)) can 
be found in Refs. [30, 31]. 

From Eq. (2.17) the cross section of double diffractive production is equal to 

a^d = ( 2 T (K; 6 )) |. (2.18) 

In Eq. (2.18) we integrated over all possible masses without any restriction. We do not expect the 
cross section for double diffractive production of small masses to be large, and we believe that most of the 
contribution in this region of masses stem from the Good Walker mechanism [49], which we consider in the 
next section. 

3 Main formulae and its phenomenological parameters 
3.1 Two channel approximation 

In the previous section we reviewed the theoretical input from the CGC/saturation approach, used for 
calculating the Green function of the resulting Pomeron. In this section we discuss a model approach 
making two simplifications: one, that the eikonal formula for the hadron scattering amplitude satisfies 
s-channel unitarity; and two, that the simplified two channel model describes diffractive production in the 
low mass region. In this model, we replace the rich structure of the diffractively produced states, by a single 
state, with wave function The observed physical hadronic and diffractive states can then be written 

iph = a'ki + ,0 4/2 ; V’n = —/3 Ti-|-a ^ 2 ; where = 1 . (3.1) 

Functions 'ipi and i /'2 form a complete set of orthogonal functions {tpi} which diagonalize the interaction 
matrix T 

-^i,k . (3-2) 
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The unitarity constraints can be written as 


(s, 6 ) 


\-^i,k ('S) b) 1 ^ + G'^i}c{s, b), 


(3.3) 


where G"^. denote the contributions of all non diffractive inelastic processes, i.e. it is the summed probability 
for these final states to be produced in the scattering of a state i off a state k. In Eq. (3.3) ^/s = W, is 
the energy of the colliding hadrons, and b denotes the impact parameter. A simple solution to Eq. (3.3) at 
high energies, has the eikonal form with an arbitrary opacity rijfc, where the real part of the amplitude is 
much smaller than the imaginary part. 

b) = i(^l- exp ^ ^ ; (3.4) 


G7!kis, b) = l- exp {-Qi^kis, b)). (3.5) 

Eq. (3.5) implies, that the probability that the initial projectiles (i, k) will reach the final state interaction 
unchanged, regardless, of the initial state re-scatterings, is given by = exp (— 6 )). 

The physical observables in this model can be written as follows 


elastic 

amplitude : 

^el ('®5 ) 

i -|- 2a^/3^ Ai^2 + /3^^2,2) ; 

(3.6) 

elastic observables : 

O'tot ~ J 



optical 

theorem : 

2lmAi^kis 


(3.7) 

single 

diffraction : 

GW 1 

-J 

' dH {aP{-a^Ai,, + {a^ - + ld‘^A2,2}f ; 

(3.8) 

double 

diffraction : 

GW 

^dd — 

f d?b a^fi^{Ai^i — 2Ai^2 + ^ 2 , 2 }^- 

(3.9) 


We denote by ‘GW’ the Good -Walker component, that is responsible for diffraction in the small mass 
region. 

In the eikonal approach we parametrize the arbitrary functions fljfc (s, b) in the form 


^ik{s,b) = J (fb'(Pb"gi{mi,b') gk{mkb") 

where the vertex gi {rrii, b') is parameterized as 




(r(y, 6 - 6 '-r)), (3.10) 


gi[mi,b') = giSip {mi,b') , where Sjp {mi,b') = —mfb'Ki{mib') 


(3.11) 


S]p ( 6 ) in Eq. (3.11) is the Fourier transform of the form factor 1/(1 -|- and Ki (z) denotes the 

modified Bessel function of the second kind (the McDonald function, see formulae 8.4 in Ref. [29]). 
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Figure 3. The net diagram that gives the largest contribution if gi ^ Gsp. The wavy lines describe BFKL 
Pomerons. The blobs denote the triple Pomeron vertices. 

3.2 Small phenomenological parameters and net diagrams 

The fits to experimental data (see Refs. [6, 16] and references therein) led to an unexpected result viz: the 
value of vertices for Pomeron-hadron interactions (gi), extracted from the fit of the data, turn out to be 
much larger than the vertex of triple Pomeron interaction ( see Fig. 3). 

Due to this small parameter G^jp/gi ^ 1, we can improve the eikonal approximation, and sum a more 
general class of diagrams, the so called net diagrams shown in Fig. 3. The analytical expression for the sum 
of the ‘net’ diagrams has been found (see Ref. [16]). Using the inequality observed in our previous paper[6]: 
viz. m in Eq. (2.10) turns out to be much larger than both mi and m 2 - Assuming this inequality to hold, 
we can simplify the integrals over h' and h” in Eq. (3.10), re-writing this equation in the form 

^ik (s, b) = J dH'd^b”gi {mi, h') g^ {mkh”) G" (t (t, h-h' -h'')) 

= (^j d%'' G“ (T (y, 6"))^ j d^b'g, {mu b') g^ (m^, 6 - ?') . (3.12) 


The summation of the ‘net’ diagram, is then given by the following simplified expression (see Ref. [16] 
for details) 


D (U; b) = [ d^b' - 

J I + 1. 


9 i{b') gk{b-b') G“(r) 


-b (T) 


9i{b'] +gk[b-b' 


where 




(T) = j 


= / d^bG^'^^^^^'^ {T {Y,b)). 


(3.13) 

(3.14) 

(3.15) 


The coefficient 1.29 results from the extraction of the value of Gsjp from the CGC/saturation approach, 
that has been considered in the previous section. 










3.3 Diffraction production in the region of large mass 


In this section we also include, in the process of diffraction production, the mechanism of production that 
originates from the dressed Pomeron, and has been discussed in section 2. 

For single diffraction the large mass contribution can be written 


^largemass _ ^ 

+_g 2/3^a^ 2 + /3® ^|f2,2 ■ 


where 


^f!kiy;b) = 


■ b) = [ dH' — 

■’ (1 


9i(b') gk(b-b'] G 


/ \ /^dressed 


(T) 


+ (T) 


gi(b'] +gk{b-b' 


(3.17) 


Al^^i{Y,Ymax,Y^in]b) = (3.18) 

/dV «„(>', Sr{T,-b\,n,) Sr{T,-V,,n) . 

where adiff{Y,Ymax,Ymin,r) is given by Eq. (2.14), Y = ln(s/so) , Ymax = In /sq) and Ymin = 

In (lVfmm/'®o)- Admax and Mmin are the largest and smallest mass produced in the diffractive processes. 

Eq. (3.16) has simple physical meaning: each term is the product of probability to produce a mass 
diffractively from the dressed Pomeron (term exp(—^fl)), and the probability of the process of single 
diffraction, from the dressed Pomeron [Ai-k^i). 

For the double diffraction production at large mass we have 

= j + 2 a^_l_ g-2nn2(>";f') I ^ (3.19) 

Afj, = j d^bgigkS^^j^{b)(Jdd{Y) where 5)^^ (6) = J d?b'Sp[b',mi) Sp{b-b\m^ , (3.20) 

where add is given by Eq. (2.18). 


3.4 Estimates for the values of the phenomenological parameters 

We have two sets of phenomenological parameters, which need to be determined by fitting to the experi¬ 
mental data. The first set is related to the description of the dressed Pomeron: ^q, A and m. The parameter 
(pQ, in principle, could be found from the solution of the BFKL equation at low energy. Unfortunately, we 
do not know the initial condition for BFKL evolution. This is the reason that we extract this parameter 
by fitting to the experimental data. The value of 4>q should be of the order of as, and therefore, we expect 
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model 

A 

4>o 

91 {GeV-^) 

92 {GeV-^) 

m{GeV) 

mi{GeV) 

777-2 (Ge 17) 

/3 

2 channel 

0.38 

0.0019 

110.2 

11.2 

5.25 

0.92 

1.9 

0.58 

1 channel 

0.323 

0.019 

25.7 

n/a 

6.35 

0.813 

n/a 

n/a 


Table 1. Fitted parameters of the model. 


that this parameter will be small. A determines the energy dependence of the saturation scale. We know 
the theoretical value of A = 26is(^ (1) — '0 (|) ^/(l ~ 7cr) where ipiz) denotes the Euler psi-function (see 
formulae 8.36 of Ref. [29]) and 7cr- = 0.37 (see Ref.[l]). On the other hand, the value of A ~ 0.3 has been 
extracted from DIS and nucleus-nucleus scattering for the energy dependence of the saturation scale[33-35]. 
The mass m plays a twofold role: it determines the impact parameter dependence, and it gives the size of the 
typical dipole in a hadron. This is a dimensional parameter, whose value is determined by non-perturbative 
QCD, and, at the moment, we have no theoretical input for this quantity. In all our formulae, we use the 
intuitive assumption that the mass m, is the largest mass in our model. 

The second set of parameters: gi and m* as well as values of a, are associated with the description of the 
wave functions in our two channel model. This set is of non-perturbative origin, and has to be determined 
by htting to the data. In our formulae we assume that ^ and m ^ rrii. This assumption is based 
on our past experience with soft Pomeron models for high energy scattering [16]. 


4 The result of the fit 
4.1 Cross sections 

We have eight parameters to be determined by fitting to the experimental data on total, inelastic and 
elastic cross sections, single and double diffractive production cross sections, and the slope of the forward 
elastic differential cross section. The value of the minimal energy for data that we use is VP = 0.546 TeV, as 
starting from this energy the CGC/saturation approach, is able to describe the data on inclusive production 
in proton-proton collisions (see Ref. [33]) . For lower energies, saturation occurs in ion-ion and proton-ion 
collisions, but not in proton-proton collisions [34]. 

The quality of the fit can be judged from Fig. 4. The fitted parameters are tabulated in Table I. 

From Table 1 we see that the qualitative features of the two channel model are similar to those of 
the one channel model [6]. The values of all the parameters have the same hierarchy i.e. g^ ^ G^jp ( 
G^jp = 1.29 GeV~^ in our approach) and m S> mj. We recall that this hierarchy is necessary for our 
formulae to be valid. A ~ 0.3, follows from the estimates of the energy dependence of the saturation 
momentum from DIS and ion-ion scattering experiments. (pQ is small, as is expected from QCD. 

We obtain a good description of the data, including the energy behaviour of the elastic slope and the 
single diffraction cross section. Therefore, we have elucidated the problems that occurred in the one channel 
model[6]. We show in this paper, that the oscillatory behavior of the single diffraction, as well as the weak 
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energy dependence of the elastic slope were artifacts of the oversimplified one channel model, and are not 
associated with our main theoretical input: the CGC/saturation approach. 

A deficiency in our two channel model is the description of the double diffractive cross section, as can be 
seen in Fig. 4. This failure is due to the small Good-Walker contribution for the production of low masses. 
It predicts a cross section of about 1 mb. To understand the sources of this small value, we note that 
Eq. (3.10) leads to the cross section for double diffraction, that can be calculated using the factorization 
relation (see Ref. [14] for example). 

aw 

Using our parameters in Table I we found that Bgd ~ B^^i ~ Bdd, where B are the slopes of the various cross 
sections. Based on our estimates for elastic and single diffraction cross section, we obtain ~ 3m6 at W 
= 7 TeV. If we evaluate the cross section of Good-Walker component using Eq. (4.1) we obtain reasonable 
description of the experimental data (see dotted curve in Fig. 4-f). We wish to stress that the shadowing 
corrections that follows from the CGC/saturation approach, do not violate the factorization properties of 
Eq. (4.1). The source of the violation of Eq. (4.1), is the net diagrams, which in our model incorporate 
the interactions of the dressed Pomerons. To illustrate the influence of this interaction, we calculated the 
slopes at W = 7 TeV. We obtain : B^i = 19.45 GeU“^; Rsd = 25.08 GeU“^ and Bdd = 44GeU“^. Using 
Eq. (4.1) with these values, decreases the estimated add by a factor of almost three . Hence, the main effect 
of the dressed Pomerons interactions is to increase the typical 6, in the impact parameter distribution for 
double diffractive cross sections. 


On the other hand, the small cross section for double diffraction arises naturally in other theoretical 
approaches [14], and the cross section, measured by the TOTEM collaboration [40] in the restricted region 
of the produced mass is rather small {add = 116 ± 25 fib for AYm = 1-8 and \Y — Y^axl = 4.7). Bearing 
this in mind, we are not too disappointed with our description of the double diffractive cross section. 

In Table 2 we present the results of our model for the various cross sections at different energies. Note, 
that for single and double diffraction, the large mass region is responsible for substantial contributions (more 
than half of the cross section). These contributions originate from the structure of the dressed Pomeron, 
and are the main theoretical input from the GGG/saturation approach. 


4.2 Amplitudes 

It is instructive to compare the two models which are both based on GGG/saturation approach: the one 
channel model developed in Ref. [6] and the two channel model of this paper. In Fig. 6 we plot the amplitudes 
Aik the two channel model and elastic amplitude for the one channel model, at different energies, as 
functions of the impact parameter, b. From Fig. 6-a we note that one of our amplitude (^ 2,2 is rather 
small at 7 TeV, while the amplitude An reaches the unitarity limit for b < 1.5 fm. It is interesting to 
compare this behaviour with the results of two channel model, based on the Pomeron as a Regge pole [16]. 
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w 

(TeV) 

^tot 

(mb) 

cjez(mb) 

(mb) 

Bel 

(GeR- 2 ) 

single 

(mb) 

diffraction 

(mb) 

double 

(mb) 

diffraction 
<^dd^ (mb) 

0.576 

62.3 

12.9 

15.2 

5.64 

1.85 

0.7 

0.46 

0.9 

69.2 

15 

16 

6.25 

2.39 

0.77 

0.67 

1.8 

79.2 

18.2 

17.1 

7.1 

3.35 

0.89 

1.17 

2.74 

85.5 

20.2 

17.8 

7.6 

4.07 

0.97 

1.62 

7 

99.8 

25 

19.5 

8.7 

6.2 

1.15 

3.27 

8 

101.8 

25.7 

19.7 

8.82 

6.55 

1.17 

3.63 

13 

109.3 

28.3 

20.6 

9.36 

8.08 

1.27 

5.11 

14 

110.5 

28.7 

20.7 

9.44 

8.34 

1.27 

5.4 

57 

131.7 

36.2 

23.1 

10.85 

15.02 

1.56 

13.7 


Table 2. The values of cross sections versus energy. and denote the cross sections for diffraction 

dissociation in the low mass region, for single and double diffraction, which stem from the Good-Walker mechanism. 
While (T^J^ and are used for diffraction in high mass, coming from the dressed Pomeron contributions. 


One notes a drastic difference in the dependence on the impact parameter, as well as the relative values 
of the amplitude. Qualitatively, the Pomeron interaction leads to stronger shadowing corrections, than the 
CGC/saturation approach. 

We would like to draw the reader’s attention to the energy behavior of the ‘black’ Ai^i amplitude, as 
shown in Fig. 5. We note that the size of the black part of the amplitude for our model, is much smaller than 
for the Pomeron based model, discussed in [16]. The second observation is that this amplitude becomes 
transparent at low energies . One of the unpleasant features of the Pomeron (Regge pole) models is the 
fact that the black component already appears at very low energies. The present model does not have this 
deficiency, providing a smooth transition from transparent to the black disc picture. 

Comparing Fig. 6 -c and Fig. 6 -d we conclude that the two channel model leads to weaker shadowing 
corrections than the one channel model. 

Looking at Fig. 6 , it is difficult to avoid the pessimistic conclusion, that until we know the non pertur¬ 
bative QCD structure of the hadrons, the high energy interaction of hadrons cannot be treated in a unique 
fashion, but is doomed to be sensitive to ad hoc assumptions. 

This statement can be illustrated by Fig. 7 where we plot the da/db"^ for the diffraction production 
processes. The low mass diffraction has pronounced peripheral features, having minima at b=0 and at b = 
1 1.5 fm. However, for high mass diffraction the main contribution stems from b=0. Single diffraction 

in the region of high mass does not have a minimum at b= 0 , while in double diffraction such a minimum 
is seen, but it is very shallow. Comparing these distributions with the one channel model, we see that 
for single diffraction, the pattern of the impact parameter distribution is similar in both cases. Double 
diffraction has a very clear peripheral distribution in the one channel model, but in the two channel model, 
the double diffraction dependence on b is similar to the single diffraction. We believe that this comparison 
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shows, that double diffraction is sensitive to the production mechanism and, in particular, to the model of 
the hadron structure. 

We would like to point out the similarity of the distributions of our two channel dipole model (shown 
in Fig.7), with those obtained by the Durham group [14], inspite, of the difference in the basic premises of 
the two groups. 

5 Conclusions 

In this paper we demonstrate that the generalization of the CGC/saturation based model to two chan¬ 
nels, allows us to describe the experimental data, without having an oscillating behavior for the single 
diffraction production, and predicting a value for the elastic slope at W—7 TeV, which is smaller than the 
experimentally measured value, ft should be stressed, that the values of the phenomenological parameters 
extracted from the fit satisfy the theoretical expectations. In addition we find new phenomenological small 
parameters: viz. G^wldi ^ 1 and mijm ^ 1 , which allows us to simplify the theoretical formulae. 

As has been discussed above, our approximation is valid only for energy less than (see Eq. (2.1)) 

Y = ln(s/so) < Vmax = — In —) « 29 (5.1) 

^BFKL ^ ^BFKL ^ 

for A = 0.38 (see Table 1). Therefore, we can use our approach for energies up to VF = 100 TeE. 

We believe that this paper together with Ref. [6], lends support to the supposition, that a consistent 
model, based on the BFKL Pomeron and the CGC/saturation approach, can be built. We have demon¬ 
strated that this model successfully describes data for high energy hadron scattering. In addition, we hope 
that this paper provides credence to the arguments, that the matching with long distance physics, (where 
the confinement of quarks and gluons is essential), can be reached within the CGC/saturation approach, 
this, without requiring that the soft Pomeron should appear (as a Regge pole). 

6 Acknowledgements 

We thank our colleagues at Tel Aviv university and UTFSM for encouraging discussions. Our special thanks 
go to Carlos Contreras, Alex Kovner and Misha Lublinsky for elucidating discussions on the subject of this 
paper. This research was supported by the BSF grant 2012124 and the Fondecyt (Chile) grant 1140842. 

References 

[1] Yuri V Kovchegov and Eugene Levin, “ Quantum Choromodynamics at High Energies ", Cambridge 
Monographs on Particle Physics, Nuclear Physics and Cosmology, Cambridge University Press, 2012 . 

[2] E. A. Kuraev, L. N. Lipatov, and F. S. Fadin, Sov. Phys. JETP 45, 199 (1977); Ya. Ya. Balitsky and L. N. 
Lipatov, Sov. J. Nucl. Phys. 28, 22 (1978). 


13 




[3] L. N. Lipatov, Phys. Rep. 286 (1997) 131; Sov. Phys. JETP 63 (1986) 904 and references therein. 

[4] P.D.B. Collins, "An introduction to Regge theory and high energy physics", Cambridge University Press 1977. 

[5] J. Bartels, C. Contreras and G. P. Vacca, arXiv: 1411.6670 [hep-th]. 

[6] E. Gotsman, E. Levin and U. Maor, Eur. Phys. J. C 75 (2015) 1, 18 [arXiv:1408.3811 [hep-ph]]. 

[7] A. Donnachie and P.V. Landshoff, Nucl. Phys. B231, (1984) 189; Phys. Lett. B296, (1992) 227; Zeit. Phys. 
C61, (1994) 139. 

[8] E. Gotsman, E. Levin and U. Maor, Eur. Phys. J. C 71, 1553 (2011) [arXiv:1010.5323 [hep-ph]]. 

[9] E. Gotsman, E. Levin, U. Maor and J. S. Miller, Eur. Phys. J. C 57, 689 (2008) ]arXiv:0805.2799 [hep-ph]]. 

[10] A. B. Kaidalov and M. G. Poghosyan, arXiv:0909.5156 ]hep-ph]. 

[11] A. D. Martin, M. G. Ryskin and V. A. Khoze, arXiv:1110.1973 ]hep-ph]. 

[12] S. Ostapchenko, Phys.Rev. D 81, 11402 (2010). 

[13] E. Gotsman, E. Levin and U. Maor, Phys. Lett. B 716 (2012) 425 [arXiv: 1208.0898 ]hep-ph]]. 

[14] V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 74 (2014) 2756 ]arXiv:1312.3851 ]hep-ph]]. 

[15] Phys. Rev. D 83, 014018 (2011) ]arXiv:1010.1869 ]hep-ph]]. 

[16] E. Gotsman, E. Levin and U. Maor, “A comprehensive model of soft interactions in the LHC era,” 
arXiv:1403.4531 ]hep-ph]. 

[17] Y. V. Kovchegov and E. Levin, Nucl. Phys. B 577 (2000) 221 ]hep-ph/9911523]. 

[18] A. H. Mueller and B. Patel, Nucl. Phys. B425 (1994) 471. A. H. Mueller and G. P. Salam, Nucl. Phys. B475, 
(1996) 293. ]arXiv:hep-ph/9605302]. G. P. Salam, Nucl. Phys. B461 (1996) 512; E. lancu and A. H. Mueller, 
Nucl. Phys. A730 (2004) 460 ]arXiv:hep-ph/0308315]; 494 ]arXiv:hep-ph/0309276]. 

[19] T. Altinoluk, C. Contreras, A. Kovner, E. Levin, M. Lublinsky and A. Shulkim, Int. J. Mod. Phys. Conf. Ser. 
25 (2014) 1460025; T. Altinoluk, N. Armesto, A. Kovner, E. Levin and M. Lublinsky, JHEP 1408 (2014) 
007; T. Altinoluk, A. Kovner, E. Levin and M. Lublinsky, JHEP 1404 (2014) 075 ]arXiv:1401.7431 
]hep-ph]].; T. Altinoluk, C. Contreras, A. Kovner, E. Levin, M. Lublinsky and A. Shulkin, JHEP 1309 (2013) 
115. 

[20] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rep. 100 (1983) 1. 

[21] A. H. Mueller and A. 1. Shoshi, Nucl. Phys. B 692 (2004) 175 ]arXiv:hep-ph/0402193]. 

[22] 1. Balitsky, ]arXiv:hep-ph/9509348]; Phys. Rev. D60, 014020 (1999) ]arXiv:hep-ph/9812311]; 

Y. V. Kovchegov, Phys. Rev. D60, 034008 (1999), ]arXiv:hep-ph/9901281]. 

[23] A. H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427. 

[24] L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 2233, 3352; D50 (1994) 2225; D53 (1996) 458; 
D59 (1999) 09400. 

[25] M. Froissart, Phys. Rev. 123 (1961) 1053; 

A. Martin, “Scattering Theory: Unitarity, Analitysity and Crossing." Lecture Notes in Physics, 

Springer-Verlag, Berlin-Heidelberg-New-York, 1969. 

[26] A. H. Mueller and D. N. Triantafyllopoulos, Nucl. Phys. B 640 (2002) 331 
]hep-ph/0205167]; D. N. Triantafyllopoulos, Nucl. Phys. B 648 (2003) 293 ]hep-ph/0209121]. 


14 


[27] E. lancu, K. Itakura and L. McLerran, Nucl. Phys. A 708 (2002) 327 [hep-ph/0203137]; 

[28] E. Levin, JHEP 1311 (2013) 039 [arXiv:1308.5052 [hep-ph[[. 

[29] I. Gradstein and I. Ryzhik, Table of Integrals, Series, and Products, Fifth Edition, Academic Press, London, 
1994. 

[30] K. G. Boreskov, A. B. Kaidalov, V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. G 44, 523 (2005) 
[arXiv:hep-ph/0506211[. 

[31] A. Kormilitzin, E. Levin and A. Prygarin, Nucl. Phys. A 813, 1 (2008) [arXiv:0807.3413 [hep-ph[[; E. Levin 
and A. Prygarin, Eur. Phys. J. G 53 (2008) 385 [hep-ph/0701178[. 

[32] M. L. Good and W. D. Walker, Phys. Rev. 120 (1960) 1857. 

[33] E. Levin and A. H. Rezaeian, Phys. Rev. D 83 (2011) 114001 [arXiv:1102.2385 [hep-ph[[, AIP Gonf. Proc. 
1350 (2011) 243 [arXiv: 1011.3591 [hep-ph[[, Phys. Rev. D 82 (2010) 054003 [arXiv: 1007.2430 [hep-ph[[, Phys. 
Rev. D 82 (2010) 014022 [arXiv:1005.0631 [hep-ph[[; 

[34] D. Kharzeev and E. Levin, Phys. Lett. B 523, 79 (2001) [nucl-th/0108006]. 

[35] J. L. Albacete, N. Armesto, J. G. Milhano, P. Quiroga-Arias and G. A. Salgado, Eur. Phys. J. G 71 (2011) 
1705 ]arXiv:1012.4408 ]hep-ph]]. 

[36] M. G. Poghosyan, J. Phys. G G 38, 124044 (2011) ]arXiv:1109.4510 [hep-ex]]. ALIGE Gollaboration, “First 
proton-proton collisions at the LHC as observed with the ALICE detector: measurement of the charged particle 
pseudorapidity density at i/s = 900 GeV,” arXiv:0911.5430 [hep-ex]. 

[37] G. Aad et al. [ATLAS Gollaboration], Nature Gommun. 2, 463 (2011) ]arXiv:1104.0326 [hep-ex]]. 

[38] GMS Physics Analysis Summary: “Measurement of the inelastic pp cross section at aLZs = 7 TeV with the 
GMS detector", 2011/08/27. 

[39] F. Ferro [TOTEM Gollaboration], AIP Gonf. Proc. 1350, 172 (2011) ; G. Antchev et al. [TOTEM 
Gollaboration], Europhys. Lett. 96, 21002 (2011), 95, 41001 (2011) ]arXiv:1110.1385 [hep-ex]]. 

[40] G. Antchev et al. [TOTEM Gollaboration], Phys. Rev. Lett. Ill (2013) 26, 262001 ]arXiv:1308.6722 [hep-ex]]. 

[41] E. Gotsman, E. Levin and U. Maor, “A Soft Interaction Model at Ultra High Energies: Amplitudes, Cross 
Sections and Survival Probabilities,” arXiv:0708.1506 ]hep-ph]. 

[42] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Phys. Rev. D59, 014014 (1999), 
]arXiv:hep-ph/9706377]; Nucl. P/iys. B504, 415 (1997), ]arXiv:hep-ph/9701284]; J. Jalilian-Marian, 

A. Kovner and H. Weigert, Phys. Rev. D59, 014015 (1999), ]arXiv:hep-ph/9709432]; A. Kovner, 

J. G. Milhano and H. Weigert, Phys. Rev. D62, 114005 (2000), ]arXiv:hep-ph/0004014]; E. lancu, 

A. Leonidov and L. D. McLerran, Phys. Lett. B510, 133 (2001); ]arXiv:hep-ph/0102009]; Nucl. Phys. A692, 
583 (2001), ]arXiv:hep-ph/0011241]; E. Ferreiro, E. lancu, A. Leonidov and L. McLerran, Nucl. Phys. A703, 
489 (2002), ]arXiv:hep-ph/0109115]; H. Weigert, Nucl Phys. A703, 823 (2002), ]arXiv:hep-ph/0004044]. 

[43] E. Levin and A. H. Rezaeian, Phys. Rev. D 83 (2011) 114001 ]arXiv:1102.2385 ]hep-ph]], AIP Gonf. Proc. 
1350 (2011) 243 [arXiv: 1011.3591 ]hep-ph]], Phys. Rev. D 82 (2010) 054003 [arXiv: 1007.2430 ]hep-ph]], Phys. 
Rev. D 82 (2010) 014022 ]arXiv:1005.0631 ]hep-ph]]; P. Tribedy and R. Venugopalan, Nucl. Phys. A 850 
(2011) 136 [Erratum-ibid. A 859 (2011) 185] [arXiv: 1011.1895 ]hep-ph]]. 

[44] E. Levin and A. H. Rezaeian, Phys. Rev. D 84 (2011) 034031 ]arXiv:1105.3275 ]hep-ph]]; K. Dusling and 
R. Venugopalan, Phys. Rev. Lett. 108, 262001 (2012) [arXiv: 1201.2658 ]hep-ph]], “Evidence for BFKL and 


15 


saturation dynamics from di-hadron spectra at the LHC,” arXiv: 1210.3890 [hep-ph]. A. Kovner and 
M. Lublinsky , “Angular and long range rapidity correlations in particle production at high energy,” 
arXiv: 1211.1928 [hep-ph] and references therein. 

[45] E. Levin and M. Lublinsky, Nucl. Phys. A 763 (2005) 172 [arXiv:hep-ph/0501173[; Phys. Lett. B 607 (2005) 
131 [arXiv:hep-ph/0411121[; Nucl. Phys. A 730 (2004) 191 [arXiv:hep-ph/0308279[. 

[46] M. A. Braun, Phys. Lett. B632 (2006) 297 [arXiv:hep-ph/0512057[; Eur. Phys. J. C16 (2000) 337 
[arXiv:hep-ph/0001268[; Phys. Lett. B483 (2000) 115 [arXiv:hep-ph/0003004[; Eur. Phys. J. C33 (2004) 113 
[arXiv:hep-ph/0309293[; C6, 321 (1999) [arXiv:hep-ph/9706373[. M. A. Braun and G. P. Vacca, Eur. Phys. J. 
C6 (1999) 147 [arXiv:hep-ph/9711486[. 

[47] E. Levin and K. Tuchin, Nucl. Phys. B 573 (2000) 833 [arXiv:hep-ph/9908317[. 

[48] E. Levin, J. Miller and A. Prygarin, Nucl. Phys. A806 (2008) 245, [arXiv:0706.2944 [hep-ph[[. 

[49] K. J. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59 (1998) 014017 [hep-ph/9807513[. 

[50] S. Nussinov, “Is the Froissart bound relevant for the total pp cross section at s = (14TeP)^?,” arXiv:0805.1540 
[hep-ph]. 

[51] A. Kovner and U. A. Wiedemann, Phys. Rev. D 66, 051502, 034031 (2002) [hep-ph/0112140,hep-ph/0204277[; 
Phys. Lett. B 551, 311 (2003) [hep-ph/0207335[. 

[52] E. Gotsman, E. Levin and U. Maor, Eur. Phys. J. G 71, 1553 (2011) [arXiv:1010.5323 [hep-ph[[. 

[53] J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012) 

[54] K. A. Goulianos [GDF Gollaboration], Acta Phys. Polon. B 33 (2002) 3467 [hep-ph/0205217] and references 
therein. 

[55] V. A. Khoze, A. D. Martin and M. G. Ryskin, arXiv:1402.2778 [hep-ph]. 

[56] A. B. Kaidalov and M. G. Poghosyan, “Description of soft diffraction in the framework of reggeon calculus: 
Predictions for LHC,” arXiv:0909.5156 [hep-ph[. 

[57] R. Giesielski and K. Goulianos, PoS IGHEP 2012 (2013) 301 [arXiv:1205.1446 [hep-ph[[. 

[58] F. Oljemark, 15th Int.Gonf. on Elastic and Diffractive Scattering (Saariselka) Finland, September 2013; G. 
Antchev et al. (TOTEM Gollaboration, Phys.Rev.Lett.Ill, 262001, (2013). 

[59] K. Goulianos, 15th Int.Gonf.on Elastic and Diffractive Scattering (Saariselka) Finland, September 2013. 

[60] B. Abelev et al. [ALIGE Gollaboration], Eur. Phys. J. G 73 (2013) 2456 [arXiv: 1208.4968 [hep-ex]]. 


16 








Figure 4. Comparison with the experimental data: the energy behaviour of the total (Fig. 4-a), inelastic (Fig. 4-b), 
elastic cross sections (Fig. 4-c), as well as the elastic slope (i?eZ,Fig. 4-d) aird single diffractioir (Fig. 4-e) aird double 
diffraction (Fig. 4-f) cross sections. The solid lines show our present fit. The data has been taken from Ref.[53] for 
energies less than the LHC energy. At the LHC energy for total and elastic cross section we use TOTEM data[39] 
aird for single and double diffraction cross sections are taken from Ref. [36]. The dotted line in (Fig. 4-f) is based on 
Eq.(4.1), (see text following this equation). 
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Figure 5. The impact parameter dependence of the ‘black’ Ai i amplitude, as function of energy W. 


18 





b in fm 

Fig. 6-a 



Fig. 6-b 



b in fm 



Fig. 6-c 


Fig. 6-d 


Figure 6. Amplitudes for one and two channel model versus b at different energies. Fig. 6-a shows the amplitude 
Aik{b) at W=7 TeV as a function of the impact parameter b in our model. The same amplitudes appear in Fig. 6-b 
for the two channel model based on the Pomeron interaction [16]. The energy behavior of the elastic amplitude is 
plotted in Fig. 6-a and Fig. 6-c for the one and two channels models, respectively. The figure for the one channel 
model is taken from Ref. [6]. 
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Fig. 7-a: high mass single diffraction 


Fig. 7-b: low mass single diffraction 




Fig. 7-c: high mass double diffraction 


Fig. 7-d: low mass double diffraction 




Fig. 7-e: single diffraction (1 ch. model) 


Fig. 7-f: double diffraction (Ich. model) 


Figure 7. da/dlP' versus b for different processes of diffraction production and different energies. ‘hm’(‘lm’) and 
‘sd’ (‘dd’) denote diffraction in high mass (low mass) region and single (double) diffraction respectively. Fig. 7-e and 
Fig. 7-f are taken from our paper in one channel model based on CGC/saturation approach[6]. 
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